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Some binding-energy-related quantities serving as effective order parameters have been used to analyze the 
shape phase transition in the odd Sm nuclei. It is found that the signals of phase transition in the odd Sm nuclei 
are greatly enhanced in contrast to the even Sm nuclei. A further analysis shows that the transitional behaviors 
related to pairing in the Sm nuclei can be well described by the mean field plus pairing interaction model, with 
a monotonic decrease in the pairing strength G. 
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I. Introduction 

Quantum phase transitions in nuclei have attracted a lot 
of attention from both experimental and theoretical perspec¬ 
tives (HE! , since they provide new insights into understand¬ 
ing the evolution of nuclear properties. The quantum phase 
transition is not of the usual thermodynamic type, but related 
to the equilibrium shape changes in the ground state of nuclei 
at zero temperature. It is thus also referred to as the shape 
phase transition or ground state phase transition, though the 
concept can also be applied to excited states. Evidences of the 
shape phase transition in nuclei are signaled experimentally 
through a sudden change in the properties of the ground state. 
An excellent example is provided by a set of even Sm isotopes, 
of which the evolution of its properties can be identified as the 
first-order shape phase phase transition experimentally 0 
On the other hand, odd-A nuclei can be approximately con¬ 
sidered as systems with an even-even core coupled to a sin¬ 
gle valence nucleon. It is thus expected that the properties of 
odd-A nuclei should be definitely affected by the shape phase 
transition emerging along the related odd isotope or isotone 
chains. 

Currently, analyses of the phase transitions in nuclei are 
mostly focused on even-even systems and have been carried 
out in the frame of phenomenological geometrical models of 
nuclear potential yiy], or algebraic models of nuclear struc¬ 
ture 112511 since the phase transitions in the intermediate and 
heavy mass region, such as the mass number A ~ 150 re¬ 
gion, are often out of reach of the microscopic shell models. 
In addition, phase transitions in odd-A nuclei may be more 
difficult to describe due to much more complicate dynamical 
situations 112611 in contrast to the adjacent even-even species. 
However, if one only wants to emphasize some special aspects 
of nuclear phase transitions, the shell model (under some ap¬ 
proximations) is still applicable to give a microscopic analysis 
of the phase transition in the A ~ 150 region. The purpose of 
this work is to give a microscopic analysis of the shape phase 
transition in the odd Sm nuclei in terms of the effective order 
parameters and odd-even effects. 

II. Effective order parameter 


One way of addressing quantum phase transitions is to re¬ 
sort to the potential energy approach. To define phase transi¬ 
tions in theory, it is convenient to consider a schematic “Lan¬ 
dau” potential 0 written as 


V (p)=^+x[(i —pY-.v/n £>o 


(1) 


with two control parameters 0 > x > 1 and y > 0. This kind 
of potential may be formally derived from the interacting bo¬ 
son model B25fl , which has been widely used to study quan¬ 
tum phase transitions in nuclei. It can be proven that the sys¬ 
tem has a second order phase transition at x = x c = 1/2 when 
y = 0 because the minima of V (j3), V m in, and are contin¬ 
uous, but 3 is discontinuous. More generally, the system 
will show a first order phase transition as a function of x for 
any fixed value of y > o. For example, one can show that 
V, n in is continuous but 3v J^ m is discontinuous at x = x c = 1/3 
when y = 2, which indicates the first order phase transition 
occurring at x c . The potential (Q3 can be considered as a sim¬ 
plified phenomenological nuclear potential surface varying as 
the function of the deformation j3, which indicates that one 
could take j3 equ iiib = fu to be the order parameter. As shown 
for the cases considered in Fig. [Q the order parameter [i, : 
changes continuously as a function of x, with the first deriva¬ 
tive being discontinuous at x c = 1 /2 if y = 0, corresponding 
to the second-order phase transition, while /),, jumps abruptly 
from 0 to 1 at x c = 1/3 when y = 2, corresponding to the first- 
order phase transition. It seems that both (f, and Vjnin, which 
may correspond to the ground state deformation and energy 
respectively, can be used to identify the shape phase transi¬ 
tion. However, for real nuclei, things are considerably more 
complicated partly because /!,, is not an obervable. In fact, 
instead of fi e , the so called effective order parameters 0 (ob¬ 
servables sensitive to shape phase transitions occurring within 
nuclei) are often used to identify nuclear shape phase tran¬ 
sitions and in some cases even determine their orders. The 
typical effective order parameters include the isomer shifts de¬ 
fined as v = c[(r 2 )o 2 - (r 2 ) 0l ] and v' = c'[(r 2 ) 2l - (r 2 ) 0l ] 0 ] 
with c and c' being the scale parameters, the B(E2) ratio 
B(E2;(L + 2) t -* L\)/B(E2\2\ —> OQ [16j, and the energy 
ratio El { /Eq-, 127], etc. Most of the effective order parameters 
are related to the quantum numbers of excited states, which 
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FIG. 1: (Color online) he order parameter j8 e as a function of x for 
y = 0 and y = 2, respectively. 



makes it particularly difficult use them to identify the phase 
transitions in odd nuclei. In contrast, the binding-energy- 
related quantities may serve as qualified effective order pa¬ 
rameters to identify the phase transitions in both even-even 
and odd nuclei m' since their values only depend on the num¬ 
ber of nucleons, and their experimental data are also relatively 
abundant. On the other hand, the number of nucleons in nu¬ 
clei is finite and the phase transitional behavior will be muted 
due to the finiteness of the system iQ. Specifically, instead of 
a discontinuity, sudden changes or flattening may be shown 
by the effective order parameters in nuclear shape phase tran¬ 
sitions mu. 


FIG. 2: Total binding energy B(Z,N) for the Sm nuclei (taken from 
01) shown as functions of the neutron number. 


nuclei, S 2n can also serve as a qualified effective order param¬ 
eter for identifying the ground phase transitions 0. Fora 
set of isotopes, the two-neutron separation energy S 2 n may be 
rewritten as a smooth contribution that is linear with the num¬ 
ber of valence neutron pairs N p , plus a contribution from the 
deformation |26] 


Sin = —A — BN P + S(2n)def, (4) 


III. Two-neutron separation energy 


For an atomic nucleus, the most basic characteristic is 
nuclear mass or binding energy. The total binding energy 
B(Z,N) for a nucleus with proton number Z and neutron num¬ 
ber N is defined as 0 

B(Z,N) = ZAfp +NM n — M(Z,N ), (2) 

where M p is the proton mass, M n is the neutron mass, and 
M(Z,N) denotes the nuclear mass. In experiments, data about 
the binding energy B(Z,N) is also abundant in contrast to 
other observables. One may expect that the shape phase tran¬ 
sition in an isotope should be reflected by the evolution of 
the binding energy, which for the Sm isotopes is shown in 
Fig.© However, the varied energy scale of the total binding 
energy B(Z,N) in an isotope chain is too large (~ 10 2 MeV) 
so that the signals of phase transition that are expected to ap¬ 
pear around N = 90 have been completely hidden behind the 
linear behavior of B(Z,N) as shown in Fig. [2] Thus one has to 
resort to other quantities related to the binding energy in or¬ 
der to identify the shape phase transition in an isotope chain. 
For even-even nuclei, the two-neutron separation energy Sin, 
which is defined as [28] 


S 2n = B{Z,N)-B(Z,N- 2), 


(3) 


can be considered as a primary and direct signature of the 
emergence of the shape phase transition [122,[24]. For odd-A 


where A and B are the parameters. To emphasize the occur¬ 
rence of the phase transition in the Sm nuclei, the experimen¬ 
tal data of Sin for both the even and odd Sm isotopes 0 are 
shown in Fig.© where the deformation contributions S(2n)<j e f 
are also shown to reveal the odd particle (single valence nu¬ 
cleus) effects on the phase transition. It is easy to know from 
Eq. (© that the deformation contributions .S'(2/;)j e r can be ob¬ 
tained from the data by subtracting a term linear with the num¬ 
ber of valence neutron pairs N p . Concretely, the results of 
5(2«)def are obtained from the data [[29] fitted with A = —19.8 
MeV and —19.4 MeV for the even and odd Sm nuclei respec¬ 
tively, and B = 0.66 MeV according to ©. As clearly seen 
from Fig. ©a), a noticeable feature is the sudden flattening 
near the neutron number N = 90 shown by S 2lI for the even 
Sm isotopes. Based on the analysis given in |22], the sudden 
flattening indicates the first order phase transition emerging 
in the corresponding isotopes. It is even more interesting to 
find that a similar or even more pronounced change appears in 
S 2n for the odd Sm isotopes near N = 90, which indicates that 
the first-order phase transition also occurs in these odd Sm 
nuclei. As shown in Fig. ©b), the phase transition occurring 
around N = 90 corresponding to N p = 4 is explicitly mani¬ 
fested in 5(2n)d e f for both the even and odd Sm isotopes. Thus 
N p = 4 {N = 90) may be considered as the critical point of 
the phase transitions in the Sm nuclei. Particularly, the phase 
transitional signal in the odd Sm nuclei seems to be greatly en¬ 
hanced by the odd particle effect. Specifically, the amplitude 
of S(2n)d e f in the odd Sm nuclei increases about 25% near the 
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FIG. 3: The two-neutron separation energies, Si n and their deformed 
part S(2n)d e f for both the even and odd Sm nuclei are shown along 
the isotopes chain. The experimental data are taken from (29il . 


FIG. 4: (a) The odd-even mass difference D for Sm nuclei with neu¬ 
tron number N = odd fitted by the pairing model in the two schemes; 
(b) the same as in (a) but for N = even. 


critical point in comparison to that in the even Sm nuclei. 

IV. Odd-even mass difference and pairing excitation 
energy 


to give a microscopic analysis of the phase transition related 
to pairing in the Sm nuclei. Specifically, the Hamiltonian is 
written as 112811 


H = Y_ £< ( a J a i + a J a l) + Y, G >i h 


j > 


( 6 ) 


Besides the two-neutron separation energy £ 2 ,, and its de¬ 
formed part S(2n)d e f, there exist some other binding-energy- 
related quantities that can serve as effective order parame¬ 
ters to identify the shape phase transitions in odd nuclei m. 
A novel transitional signal in odd-even nuclei related to the 
binding energy B(Z,N) is given as the odd-even mass differ¬ 
ence 11711 . which is defined as 


D = B(Z,N) — 


B(Z,N-1)+B(Z,N+1) 


(5) 


where e, represents the single-particle energy of the i-th Nils¬ 
son level, and b] = ajat (b, = ajat ) is the pair creation (an¬ 
nihilation) operator with i labeling the time-reversed state of 
that labeled by i. To solve such a Hamiltonian, two ap¬ 
proximation schemes are considered for the pairing interac¬ 
tion in this work. One is the nearest-orbit pairing-interaction 
model (33 211, which is a simplified version of the Gaussian- 
type pairing interactions suitable for deformed nuclei [132 1 
with the orbit-dependent pairing strengthes written as 


Meanwhile, it is also well known that the odd-even effects, 
such as those in the odd-even mass difference, are the most 
important evidences of pairing in nuclei 1281 . Therefore, 
some microscopic factors relevant to pairing in manifesting 
the phase transition can be extracted from the odd-even ef¬ 
fects. To do that, a shell model Hamiltonian including the 
deformed mean field and pairing interaction will be adopted 


G u = a e-Pte- 6 ;) 2 , (7) 

where a < 0 and /3 > 0 are the adjustable parameters. It 
is clear that the pairing strength G\j shown in (|7} is orbit- 
dependent, and the nearer the two orbits the stronger the pair¬ 
ing interaction between the two pairs. As an approximation to 
the Gaussian-type interactions given in l[7]). only the on-orbit 
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pairing interactions G„- and the nearest-orbit pairing interac¬ 
tions G„ | | or Gn -1 are considered in the nearest-orbit pairing 
model, while G,y with |i — j\ > 2 are neglected 1 30l 31], As a 
further approximations, we set G !( - = G „+1 = G. Thus, there 
is only one free-parameter G to be determined for each nu¬ 
cleus in the isotopes. Such a pairing interaction form can be 
exactly solved for all the Sm nuclei by directly diagonalizing 
the Hamiltonian in the valence nucleon space. More details 
about the solutions_of the nearest-orbit pairing model can be 
found inRef. I3(i[ 311. For convenience, we denote the exactly 
solvable nearest-orbit pairing interaction as the ENO scheme. 
Another approximation scheme is the constant pairing interac¬ 
tion & in which the paring strength is set as G (/ - = G for all 
the single-particle orbits i. However, the constant pairing in¬ 
teractions can be exactly solved only for a nucleus with very 
few valence nucleons, due to the difficulty of computation. 
To apply the constant pairing interactions form to the whole 
chain of the Sm isotopes, the well known BCS method 12811 is 
used to solve the corresponding Hamiltonian. We denote the 
constant pairing interaction solved by the BCS theory as the 
CBCS scheme. 

By using the ENO and CBCS schemes, the odd-even mass 
difference D for the Sm nuclei has been fitted by the Nils¬ 
son mean-field plus the pairing model and the resulting values 
together with the experimental data are shown in Fig. 0 In 
theory, the odd-even mass difference is given as 


D = —E g (Z,N) 


-E g (Z,N-l)-E g (Z,N+l) 

2 


( 8 ) 


where E g (Z,N) is the ground state energy solved from © for 
a nucleus with the proton number Z and the neutron number 
N. In our calculations, the single-particle energies {e,} are 
calculated from the Nilsson model with deformation param¬ 
eters taken from |33j], which were determined systematically 
from the corresponding experimental data iHl . In addition, 
it is assumed that the odd-even mass difference in the Sm 
isotopes only comes from the neutron part since the proton 
number in a chain of isotopes is a constant. As clearly seen 
from Fig. 0 a) and (b), the experimental values of the odd- 
even mass difference in Sm can be well reproduced in the 
pairing model for both schemes. Specifically, the values of 
the odd-even mass differences D are all negative for the odd 
Sm nuclei but positive for the even Sm nuclei, which indicates 
that even-even nuclei are more bounded than the odd-even nu¬ 
clei (28]. More important, the evident phase transition signals 
in experiments shown by D, of which the values reach their 
minimum or maximum around N = 90, are nicely expressed 
by those calculated from the pairing model. It is thus con¬ 
firmed that the pairing interaction is indeed a key factor in 
driving phase transitions in nuclei. 

Further, the fitted pairing strength G in the two schemes is 
shown in Fig. 0 As shown in Fig. 0a) and (b), the result¬ 
ing pairing strength G values in both schemes show a mono¬ 
tonic decrease for the even Sm isotopes as the neutron num¬ 
ber N increases, except that the variational behavior of G in 
the ENO scheme is a little smoother than that in the CBCS 


scheme. A similar situation also appears in G for the odd Sm 
nuclei as seen in Fig. 0c) and (d). It is thus confirmed that 
the phase transition behavior related to pairing in an isotope 
may be driven by the pairing interaction with a monotonic de¬ 
crease in the pairing strength as the neutron number increases. 
In addition, one may find that the energy scale of G in the 
ENO scheme is almost ten times that in the CBCS scheme, 
as seen from Fig. 0 It is not difficult to understand this by 
considering the fact that only the on-orbit and nearest-orbits 
pairing interaction are taken into account in the ENO scheme 
for a single-particle orbit i, in contrast to the CBCS scheme. It 
seems that the the transitional behavior of the odd-even mass 
differences D can be well illustrated in theory via the shell 
model Hamiltonian including only the mean-field plus pairing 
interaction. 

To further reveal the phase transitional behaviors closely re¬ 
lated to the pairing interaction, we also calculated the pairing- 
excitation energy (PEE) & which is also a quantity sensi¬ 
tive to the pairing strength G. Since the angular momentum 
projection along the third axis in the intrinsic frame is con¬ 
sidered to be a conserved quantity in the model, the pairing- 
excitation states determined by the model are thus regarded 
approximately as excited states with the same spin and parity 
as those of the ground state of a nucleus. For example, the 
PEEs in the model for even-even nuclei are considered as the 
energies of the excited 0 + state, E 0 +. In Fig. 0 the results 
of the first PEE calculated for the even Sm nuclei in the two 
schemes and the corresponding experimental data are shown. 
Notably, only the first neutron PEEs in the two schemes are 
taken to be compared with the experimental data since the pro¬ 
ton PEE may be much higher than the corresponding neutron 
PEE in the present case. In the ENO scheme, the first PEE 
may be briefly denoted as PEE = E\ — E g , where E\ and E g 
represent the energy value of the first excited state and that 
of the ground state respectively, since broken pairs have not 
been taken into account in this scheme. In contrast, the first 
PEE calculated from the CBCS scheme can be explicitly given 
as 


PEE = 2i/(e 0 — A) 2 + A 2 , 


(9) 


where A represents the Fermi energy, £o denotes the single¬ 
particle energy closest to A, and A is the so-called gap param¬ 
eter 0. All the parameters involved in (0 can be determined 
by the standard BCS theory 1 ,2811 . As clearly seen in Fig.0 the 
first PEE in experiments also provides an evident phase transi¬ 
tional signal around N = 90. More importantly, such a phase 
transitional characteristic shown by the PEE can be well re¬ 
produced by the results calculated from the pairing model in 
the ENO scheme. It is thus further confirmed that the tran¬ 
sitional characteristics in connection with pairing in the Sm 
isotopes are indeed driven by the pairing interaction with a 
monotonic decrease in the pairing strength. However, it can 
be also found that the PEE obtained from the constant pairing 
interactions in the CBCS scheme are much higher than those 
found experimentally, and the global behavior of the PEE in 
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FIG. 5: The evolution of the pairing strength G changing as a function of neutron number. 
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FIG. 6: The first pairing-excitation energy calculated by the two 
schemes compared with that determined by experiments I29 : l. 


theory is also completely different from that in experiments. 
As a consequence, the exact solutions are important to explore 
the phase transition in the pairing model. It should be noted 
that the PEEs in the odd Sm nuclei are not taken into account 
here because some single-particle excitations with spin and 
parity the same as those of the ground state are often involved 
in the low-lying spectrum, which makes it difficult to pick out 
the PEE from the spectrum of the odd Sm nuclei. 


VI. Summary 


In conclusion, we have made a microscopic analysis of the 
shape phase transition in the odd Sm nuclei from the point of 
view of the effective order parameter. Through analyzing the 
two-neutron separation energy, it is confirmed that the first or¬ 
der phase transition also occurs in the odd Sm isotopes as it 
does in the even Sm isotopes but with the signals of the phase 
transition in the odd species greatly enhanced by the odd neu¬ 
tron effect. It is also shown that the odd-even mass differ¬ 
ences may reach their extreme value around the critical point, 
thus serving as a valid effective order parameter for the iden¬ 
tification of the shape phase transition in the odd Sm nuclei. 
Particularly, analysis based on the mean-field plus pairing in¬ 
teraction Hamiltonian shows that the critical phenomena rele¬ 
vant to pairing in the Sm nuclei can be driven by the pairing 
interaction with a monotonic decrease in the pairing strength 
G. In addition, the results also indicate that the exactly solv¬ 
able models are important to analyze transition characteristics 
related to pairing in the excited states. Although the discus¬ 
sion in this work provides a specific example of microscopic 
analysis of the shape phase transition in odd-A nuclei, inves¬ 
tigations based on a more realistic shell model Hamiltonian 
with different truncation schemes |6. 35, 361 are still needed 
to eventually confirm or disprove the theoretical predictions. 



















6 


Support from the Natural Science Foundation of China 
(11375005,11005056,11175078 and 11405080). 


[1] Iachello F. Phys. Rev. Lett., 2000, 85: 3580-3583 

[2] Iachello F 2001 Phys. Rev. Lett. 87: 052502 

[3] Jolie J, Casten R F, Brentano P von, and Werner V. Phys. Rev. 
Lett., 2001, 87: 162501 

[4] Jolie J, Cejnar P, Casten R F, Heinze S, Linnemann A, and 
Werner V. Phys. Rev. Lett., 2002, 89: 182502 

[5] Iachello F and Zamfir N V. Phys. Rev. Lett., 2004, 92: 212501 

[6] Sun Y, Walker P M, Xu F R, and Liu Y X. Phys. Lett. B, 2008, 
659: 165-169 

[7] Mu L Z and Liu Y X. Chin. Phys. Lett., 2005, 22: 1354-1357 

[8] Zhang J F, Long G L, Sun Y, Zhu S J, Liu F Y, and Jia Y. Chin. 
Phys. Lett., 2003, 20: 1231-1233 

[9] Zhang D Land Liu YX. Chin. Phys. Lett., 2003, 20: 1028-1030 

[10] Zhang D L and Zhao H Y. Chin. Phys. Lett., 2002, 19: 779-781 

[11] Zhang Z, Zhang Y, An Y, Li S X, and Xu J. Chin. Phys. Lett., 
2013, 30: 102101 

[12] Dai L R, Teng W X, Pan F, and Wang S H. Chin. Phys. Lett., 
2011,28: 52101 

[13] Warner DD and Casten RF. Phys. Rev. C, 1983, 28: 1798-1806 

[14] Liu Y X, Mu L Z, and Wei H. Phys. Lett. B, 2006, 633: 49-53 

[15] Zhang Y, Hou Z F, Chen H, Wei H, and Liu Y X. Phys. Rev. C, 
2008, 78: 024314 

[16] Zhang Y, Hou Z F, and Liu Y X. Phys. Rev. C, 2007, 76: 
011305(R) 

[17] Zhang Y, Bao L N, Guan X, Pan F, and Draayer J P. Phys. Rev. 
C, 2013, 88: 064305 

[18] Zhang Y, Pan F, Liu Y X, Luo Y A, and Draayer J P. Phys. Rev. 
C, 2013, 88: 014304 


[19] Luo Y A, Pan F, Wang T, Ning P Z, and Draayer J P. Phys. Rev. 
C, 2006, 73: 044323 

[20] Li Z P, Niksic T, Vretenar D and Meng J. Phys. Rev. C, 2009, 
80 : 06130HR) 

[21] Pan F, Draayer J P, and Luo Y. Phys. Lett. B. 2003, 576: 297- 
302 

[22] Casten R F and McCutchan E A. J. Phys. G, 2007, 34: R285- 
R320 

[23] Cejnar P and Jolie J. Prog. Part. Nucl. Phys., 2009, 62: 210-256 

[24] Cejnar P, Jolie J, and Casten R F. Rev. Mod. Phys., 2010, 82 : 
2155-2212 

[25] Iachello F and Arima A 1987 The Interacting Boson Model 
(England: Cambridge University). 

[26] Iachello F, Leviatan A, and Petrellis D. Phys. Lett. B, 2011, 705 : 
379-382 

[27] Bonatsos D, McCutchan E A, Casten R F, and Casperson R J. 
Phys. Rev. Lett., 2008, 100 : 142501 

[28] Ring P and Schuck P, The Nuclear Many-Body Problem 
(Springer-Verlag Berlin 1980) 

[29] LBNL Isotopes Project Nuclear Data Dissemination Home 
Page. Retrieved June 14, 2012. <http://ie.lbl.gov/toi2003> 

[30] Pan F and Draayer J P. J. Phys. A, 2000, 33: 9095-9100 

[31] Chen Y Y and Pan F. High Energy Phys. and Nucl. Phys., 2002, 
26 : 1142-1150 

[32] Molique H and Dudek J. Phys. Rev. C, 1997, 56: 1795-1813 

[33] Moller P, Nix J R, Myers W D and Swiatecki W J. Atomic Data 
Nucl. Data Tables, 1995, 59: 185-381 

[34] Audi G, Bersillon O, Blachot J and Wapstra A H. Nucl. Phys. 
A, 1997, 624: 1-124 

[35] Zhao Y M, Yoshinaga N, Yamaji S, and Arima A. Phys. Rev. C, 
2000, 62 : 014316 

[36] Zhao Y M, Yamaji S, Yoshinaga N, and Arima A. Phys. Rev. C, 
2000, 62 : 014315 



